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Introduction : SIS epidemic model 



Infectious diseases are those caused by pathogens (such as viruses, 
bacteria, epiphytes) or parasites (such as protozoans, worms), and 
which can spread in the population. The impact of infectious 
diseases on human and animal is enormous, both in terms of 
physical and mental sufferings, and social and economic 
consequences. 
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Different mathematical models exist to study the progress of an 
epidemic in a population: 

SIR, SIRS, SIS etc. 
e S — >• Susceptible 
a I — > Infected 



Example : SIS models infections such us common flu 
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Classic SIS model 



I dS — (fx — /xS — (3SI + 7 /) dt , 
\ dl = (-(v + ^l + (3SI)dt, 

where : 

a /x — >• birth rate (coinciding with death rate) 

9 /3 infection coefficient 
9 7 — >■ recovery rate 
Constant population =>■ S + I = 1 
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Stochastic SIS model 



a Model proposed by Gray et al. (2011) 
The perturbation : /? (3 + adB(t) 



SDE 



dS = ([i — /iS — PSI + 7 /) dt — crSIdB, 
dl = (-(/i + 7 )/ + (3 SI) dt + aSIdB , 



( 2 ) 



a Model proposed in this presentation 
Stochastic perturbation+jumps: 



P P + adB(t)+ I H(z)N(dt, dz) 
J Y 
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Where : 

a N : represents a compensated Poisson measure 
a H : affects of random jumps on f5 
a Y : measurable subset of [0, oo) 
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So, equation 2 becomes : 

( d S(t) = (M - vS(t) - /3S(t)l(t) + 7 /(t)) dt - o-S(t)/(t)de(t) - H(z)S(t)l(t)N(dt, dz), 



( dl{t) = (-(/i + 'y)l(t) + 0S(t)/(t)) dt + trS(t)l(t)dB(t) + J/ H(z)S(t)l(t)l\l(dt, dz), 

What are the advantages of using jumps ? 

1- In reality, the infection coefficient /3 can change dramatically its 
value (for example in a population with common flu a surprising 
change of weather affects quickly the value the spread of 
infection). So, the model with jumps is more suitable. 2- Bao et al. 
(2011) proved that convergence rate using jumps has order 1/p for 
any p > 2 instead of 1/2 for any p > 2 without jumps. 



(3) 
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Well posed ness 



Let define the two following sets : 

R+ = {(xi,x 2 )|x/ > 0, / = 1,2} A = jx E K+; xi +x 2 = lj , 

The following theorem ensures that our model is well posed. 
Theorem 

Define C(z,x) = (1 — H(z)x) (1 + H(z)(l — x)), for 
(z,x) € Yx (0,1). If 

sup [ f log 1 ~ ')i/(dz)< oo. (4) 

0<x<iJy\ C(z,x) J 

Then, the set A is almost surely positively invariant by the system 
(3), that is, if(S( 0), 1(0)) G A, then P((S(t), l(t)) e A) = 1 for 
all t > 0. 
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Extinction of disease 



Since we have S(t) + l(t ) = 1 we will focus on l(t ) and study the 
following SDE : 



dl{t) = l(t) (/3(1 — l(t)) — (// + 7)) eft + <r/(t)(l — l(t))dB(t) 

+ f H(z)l(t){l- l(t))N(dt,dz). (5) 

J Y 

- fi(l(t))dt + f 2 (l(t))dB(t) + f f 3 (l(t),z)N(dt,dz). 

J Y 
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Theorem 



Let assumption 4 hold and assume that 

sup [ (log(l + H(z)x)) 2 v{dz) < oo. (6) 

0<x<l J Y 



and 

If B < 0. then lim sup - log(/(t)) < C a.s. 

t—> oo t 

1 B 2 

If B > 0, then lim sup - log(/(t)) < — - + C a.s., 
t->o o t 4 A 

where 

A =\(° 2 + \ S Y H 2 {z)n{dz)) , 

B = -/3 + cr 2 + \ f Y H 2 (z)n{dz), 

C = P-(p + l+\(cr 2 + §/ Y H 2 {z)n{dz)^ . 



The theorem leads to the following corollary (extinction condition) 
Corollary 

Let assumptions (4) and (6) hold. 



If (B < 0 and C < 0) or 



B 2 N 

B > 0 and — + C < 0 



then for any initial value 1(0) e (0, 1), the solution l(t ) of the 
stochastic differential equation (5) tends to zero exponentially 
almost surely. In other words, the disease dies out with probability 



one. 



Computer simulation for extinction 




Figure: Computer simulation of a single path of /(t) for the SDE model 
(5) with initial condition /( 0) = 0.4 using the parameter values : 

At = 0.014, /3 = 0.31, 7 = 0.32, cr = 0.11, H{z) = 0.1, z e K. 
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Persistence of the disease 



We consider the two following quadratic functions 1 : 

Q(x) = -Ax 2 + Bx+ Cand Q'(x) = - A'x 2 + B'x + C' 
where 

A> = \ (v 2 + j Y H2 ( z M dz )J > A 

B' = — (3 + cr 2 + I H 2 (z)is(dz) > B , 

J Y 

C' = P- (d+l+l, (cr 2 +j Y H 2 {z)v{dz)J^ < C. 
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Theorem 



Let assumptions (4) and (6) hold. 

(/) If 0 < H(z) < 1, z E Y and C > 0 then, 

lim sup /(t) > a.s. (7) 

t—tO O 

(//') If C > 0 then, 

lim inf l(t) < £, a.s., (8) 

t-s>oo v ' 

where, £ and £' are the positives roots of Q(x) = 0 and Q'(x) = 0 
respectively. 



Computer simulation for persistence 




Figure: Results of one simulation run of SDE (5) with initial condition 
1(0) = 0.4 using the parameter values : 

H= 0.017, £=1.3, 7=0.5, cr = 0.2, H(z) = 0.95, zeK 
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Conclusion 



This work studied a stochastic SIS epidemic model with constant 
population size driven by Levy noise. We first proved the existence 
and uniqueness of its global positive solution. Then, the long-term 
behavior of the stochastic SIS epidemic model is investigated. We 
obtained sufficient criteria for extinction and persistence of the 
epidemic in the population. 
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